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1 Introduction. 

Recently some meteorologists have attempted to obtain a mathematical description of those large scale 
atmospheric motions called long waves. The primary difficulty is that the equations of gas dynamics, which 
govern the motion of the atmosphere, are so complicated that they have not been satisfactorily solved 
analytically, even when viscosity, heat-conduction and moisture are neglected. They are also unsatisfactory 
for numerical solution because of the extremely short time intervals which they necessitate. 

In attempting to simplify the equations, meteorologists have often observed that omission of all acceler- 
ation terms from the equations of motion lead to the hydrostatic pressure and geostrophic wind equations. 
The first of these results is considered to be in extremely close agreement with observation, while the second 
result is in fairly close agreement for large scale motions, particularly at high altitudes where topographical 
effects are unimportant. 

Although such a derivation of these results is logically somewhat unsatifactory, a more serious difficulty 
arises from the attempt to combine them with the remaining equations (the conservation of mass and 
constancy of entropy equations). When the hydrostatic pressure and geostrophic wind are combined with 
these equations, it is found that the pressure at the ground is essentially independent of time. 

An atempt to overcome this difficulty has been made by J. Charney. By using the observed values of 
all the quantities entering the equations, he computes the magnitude of each term in the equations. He 
then retains the largest terms in each equation. In this way he finds that the acceleration terms in the 
equations of motion are small, and thus obtains the hydrostatic pressure and geostrophic wind. In the 
mass equation the largest terms are the two which constitute the horizontal divergence of the horizontal 
wind, but the vanishing of this divergence does not yield a new result since i it follows identically from 
the geostrophic wind equations. Therefore, another equation is obtained by differentiating and combining 
the horizontal equations of motion, retaining previously neglected terms, and eliminating the horizontal 
divergence by means of the mass equation. The largest terms in the resulting equation are then retained, as 
well as the complete entropy equation. In this way a semi-empirical deduction of the hydrostatic pressure 
and geostrophic wind equation is given, and a complete set of equations embodying them is obtained. To 
further simplify the equations additional assumptions, such as that the wind is independent of height, are 
made. 

Nevertheless it still seemed to us that a systematic mathematical derivation of the hydrostatic pressure 
and geostrophic wind equations, together with simplified mass and entropy equations , would be worthwhile. 
A method which has been used to derive the shallow water theory, the membrane theory of plates, and the 
theory of thin heavy jets immediately suggested itself. 

This method involves two steps. First dimensionless variables are introduced which involve a small 
parameter that stretches some coordinates and compresses others. The parameter may represent the ratio 
of s typical vertical dimension to a typical horizontal dimension of the problem. Then it is assumed that the 

'Institute for Mathematics and Mechanics Report, NYU, May 1951. Retyped June 2006 in with typos corrected, 

t Departments of Mathematics and Mechanical Engineering, Stanford University, Stanford, CA 94305-2125, USA 
tCourant Institute of Mathematical Sciences, New York University, New York, NY 10012-1185, USA 



1 



solution can be expanded as a power series in this parameter. The expansions are inserted into the equations 
and coefficients of each power of the parameter are equated to zero, yielding a sequence of equations for 
the successive terms in the solution. If an appropriate choice of dimensionless variables has been made, the 
first terms in the solution satisfy equations of the expected type. We were guided by Charney's numerical 
estimates in selecting our dimensionless variables. 

The result is a simplified set of equations for the first terms in the solution!, embodying the hydrostatic 
pressure and geostrophic wind equations. This set does not suffer from the old difficulty of yielding a time- 
independent pressure at the ground. It is also simpler than the original set, and may yield to approximate 
solutions. In fact some appromimate solutions are given in Section]^ Furthermore these simplified equations 
are more suitable for numerical solution than the original equations. 

The primary advantage of the present method of derivation is that the dervation of equations, including 
the higher order equations, is completely automatic once the change of variables has been made. Thus our 
equations are slightly different, and in fact simpler, than Charney's because the expansion scheme determines, 
for example, that a particular coefficient should be a known zero order quantity rather a sum of known and 
higher order unknown quantities. A less systematic procedure may not yield such results because the order of 
magnitude of every term is not noted. A secondary advantage of the method is that the mathematical nature 
of the approximation, and its asymptotic character, can be understood, thus leading to an interpretation of 
the accompanying boundary layer phenomena. 



2 Exact Formulation. 

We consider the motion of a non-viscous, non-heat conducting, polytropic gas around the earth. The 
equations of motion are written in Eulerian form employing spherical coordinates referred to a rotating 
cordinate system. The axis of rotation is taken to be the polar axis. The coordinates are radius r, colatitude 
6 and longitude 4>, and u,v,w are the respective velocity components. The pressure is p, the density p, and 
the angular velocity of the coordinate system is fi, which is the angular velocity of the earth. The only 
external force is that of gravity which has the components — Gi,— 6*2 and in the r, 0, directions. The 
presence of the 9 component is due to the non spherical shape of the earth and the non-symmetrical mass 
distribution. The surface of the earth is given by r = 0) and the velocity of the gas is assumed to be 
tangential to the earth at its surface. In addition p and p are assumed to approach zero as r becomes infinite. 

With these definitions, the three equations of motion and the equations of conservation of mass and 
constancy of entropy for each "particle" become: 

(1) Ut -f- uUr + vr^^ug + (r sn\9)^^wu^ = 

r-^{u^ + -u?) + 2VLw sin 9 - fl^r sin^ 9 - Gi - V 

(2) Vt + uVr + vr^^vg + {r sin 9)^^ wv^ 

= -r^^uv + r^^w'^ coi9 + 20,11] cos9 + VL^r sin9cos9 - G2 - p^^per^^ 

(3) Wt + uWr + vr^'^wg + {r sin9)^^wW(f, 

= —r^^uw — r^^vw cot 9 — 2ilusin9 ~ 2r2wcos0 — {pr sin9)~^p^ 

(4) Pt + {pu)r + {rsin9y^{pvsin9)g + {rsin9y^{pw)^ + ^ = 

r 

(5) Pt + upr -I- vr^^pg + {rsinffy^wp^ ^ ^ '"^^^ sin 6')^-^wp0] 
The boundary condition at the earth's surface is: 

(6) u - r-^vRe-{r sin 9)-'^ wR^^O, at r = i?(6l, 0). 

Equations (l)-(5) are five equations for the five functions u, v, w,p, p assuming Gi, G2 and R are known. In 
addition the initial values of the five unknown functions are assumed to be given. 
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It is convenient to introduce the "effective" components of gravity, defined by: 



(7) 51 = Gi+n^rsin^e 

(8) 52 = G2 - OVsin6lcos6'. 

The formation of the earth is such that at the surface r = R{6,(j)), the tangential component of "effec- 
tive" gravity is nearly zero, Since the earth is almost sphereical, this component is practically g2, which is 
consequently small and usually neglected in meteorology. 

3 Dimensionless Variables. 

It is convenient to introduce the new independent variable z = r — a, a denoting the mean radius of the 
earth. Then the surface of the earth is given by 2; = Z{9, cj)) = R{0, (j)) — a. 
We now introduce dimensionless variables by means of the equations: 

u = e^CoU, z = e^az, t = e^^ac^^t, Z = e^aZ, 

(9) v = t^cav, 9 = 9Q + e9, P = PoP, Co=poPo^ 
w = e'^cow, (f> — (j)o + e(j), p — pop p = flacQ^. 

Here the barred quantities are the old variables and the unbarred quantities are the corresponding dimen- 
sionless variables. The quantities 9o,(j)n,co,Po and po are constants; cq is a velocity and po,Po arc typical 
pressure and density values. The quatity e is a small dimensionless parameter which introduces a stretching 
in the scale of some quantities and a contraction in the scale of the others. The quantity may be considered 
to represent the ratio of a typical vertical dimension of the atmosphere to the radius of the earth, and it is 
therefore very small. This small parameter will later provide the basis for a series expansion of the solution. 
We also introduce the dimensionless components of "effective" gravity, Ai and A2 by the equations 

(10) ffiocg ^ = e"^Ai, g2a%'^ = e"A2. 

The factor e" makes the smallness of 52 apparent, since we assume that a is greater than 3, but is otherwise 

unspecified. 

Now, introducing equations 7-10 into equations 1-6 we have: 

(11) e^ut + e^uu^ + e^wue(l + e^z)"^ + e'^wu^il + e^z)"^(sin^)"^ 

= e'^iu^ + w'^){l + e'^z)-^ + e'^2pwsm9 - Xi - p^^P:,, 

(12) eSt -I- e^uvz + e^vve{l + e^z)"^ -|- e^wv4,{l + e^z)"^(sin^)-^ 

= -t'uvil + e^z)-'^ + t'v? cot 9 [1 + e^z)-'^ + e^2/xw; cos (9 - e"+^A2 - p" V(l -l- e^z)-^ 

(13) e^wt + e'uw, + e\tug{l + e^z)'^ + e^ww^{sm9)-^{l + e^z)"^ 

= -£'"^^^(1 + e^zy^ - e'vwcot9{l + f^zy^ - e^2p.usm9 - e^2/iwcos^ - p^(psin^)"^(l + e^z)"\ 

(14) ept + (pu), + e{pvsm9)e(sin9)-^{l + e^z)-^ + e{pw)4,{sm9)-\l + e^z)'^ + e^2pu{l + e^z)'^ = 0, 

(15) e{pp~'^)t+u{pp-y,+e{l + e^z)-Mpp'^')e + ^{pp-'')A=^^ 

(16) u = e^vZe{l+e^Z)-^ + e^wZ^{sm9)-^{l + e^Z)-^, at z = Z{e,(j)). 

4 Power Series Solution. 

To solve equations 11-15 subject to the prescribed conditions, we asume that u, v, w,p and p can be expressed 
as power series in e. Thus we assume 
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oo oo 

(17) «; = 5]eV(0,0,2,t), p = EeV(0,0,2,t), 

i=0 i=0 

oo 

Wc now insert equations 17 into equations 11-16 and equate to zero the coefficients of each power of e. 

From the coefficients of we obtain: 

(18) -pl = Aip°, 

(19) pI = 0, 

(20) pI = 0, 

(21) (p%0), = 0, 

(22) „0(pO_^^0) ^ 0^ 

(23) u° = 0, at z = 0. 

Prom equations 21 and 23 we find w° = 0, and from equations 19, 20 we have p° = p^{z,t). Thus the 
equations 18-23 are equivalent to 

(24) ^0 = 0, /=/(z,i), -p° = Aip°. 
From the coefficients of in equations 11-16 we obtain: 

(25) -pI = Xip\ 

(26) pI = 0, 

(27) pI = 

(28) pO^(p0^i)^ + pO^;0 + ^^o = 0, 

sin ^0 

(29) + = ^(p? + nV°), 

(30) =0, at 2 = 0. 
Prom the coefficients of we have: 

(31) -pI = Aip2, 

(32) pI = 0, 

(33) pI = 0, 



(34) p^^ + (p%^ + pV), + (/t;,"+p%^+pVcot^o) + ^^(pX + /'H-'oX^cot^o) = 0, 

(35) pI+uV.+uV. = ^-^{p\+u^pI+u'pI) + C^-^-^M + u'p% 

(36) = 0, at z = 0. 

Prom the coefficients of we have (from equations 11-13) 
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(37) -pI = Aip3 

(38) 2/i«;°cos6io = ^Pg 

(39) -2^.°cos^o = ^/,. 

We will not write the remaining third order equations, since they will involve additional coefficients. 
Instead we will consider the coefficients of e'^ in equations 12 and 13, which yield 

(40) v° + u^v° + v%°g + w°vl{sm Oq)-^ = 2,iw^ cos - 2^^°^ sin - ^Pg - 

(41) w;° + v}wl + v°w^e + w^w^Asm Oq)-^ 
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= -2av} sin 6*0 - 2liv^ cosOq + 2av°e siuOq - ^ +pl{p^ sin 6*0 + p°e cosOoMp^ smOo)-'^ . 

sin ^0 

5 Consequences of the Equations. 

Before atttempting to count equations and unknowns, we will simplify the equations by deducing some 
obvious consequences of them. First, by using equations 38, 39 in equation 28 we obtain 

(42) + (/«!), = 0. 
Now using equation 18 in equation 42 yields 

(43) - Ar V°t + {P^'u'). = 0. 

Integrating with respect to z and applying the boundary condition p° = at ^ = oo, we have 

(44) pO+pV=0. 

From equations 44 and 29 we find 

(45) u\ = 0. 

Using equations 30 and 45, we finally obtain 

(46) = 0. 

Then from equation 44, = 0. Thus 

(47) P°=p\z). 

Now of the 16 quantities, p^,p^,p'^,p^,p'^, p^,p^,p^, p^,u°,u^,u'^,v°,v^,'w'^,'w^, which appear in equations 
18-47, two, li" and u^, are zero (see equations 24, 46). A third, p", is independent of t (by eq. 47) and is 
therefore determined by the initial data. Of the remaining 13 quantities, 9 — , p^ , p"^ , p^ jv'^ ,w'^ 
— are given explicitly in terms of the remaining 4, p^,p^,p^,p'^. Of these 4, p'^ automatically drops out when 

and are eliminated (see eq. 48). Furthermore p"^ appears only in the equation for p^. Thus if only p^ 
and p^ can be determined, then p'^ ,p^ ,p^ , p'^ , p^ , p^ ,v'^ , w° will be known. 

To obtain equations for the determination of and p^, we first attempt to eliminate and from 
eq. 34 by means of eqs. 40 and 41. To this end we differentiate eq. 41 wih respect to (f) and divide it by 
sin^O) differentiate eq. 41 with respect to 9, and subtract the second from the first. We then obtain 



1 Iji 

(48) 2pcos9o{-r-^wl+vl) = {v^+v%^o + -: 



smOo f ' ' * ' " ' sin6'o''^sin6lo 

/ , 0, ^^e^cot^o p|cos2^o 
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The expression on the left also appears in equation 34. It is to be noted that p"^ does not occur in equation 
48. Now using eq. 48 in eq. 34, and eliminating some terms by the aid of eqs. 38, 39, we have 



0„.2\ , P r/„,0 , „,0„,0 , ^ ^ 



2/xcos6'o smfc'o sm^o 

-K° + v^w'e + - 2p|(p° cos^o)"^] = 0. 

sin (70 

Equation 49 together with equations 39, 38, 35 and 25 are five equations involving the six unknown 
fimctions v'^ , w'^ , , p"^ , . and . The only other unused equations involving any of these quantities are 
equations 26 and 27 which simply imply 

(50) p'=p\z,t). 

Thus the above equations alone do not seem adequate for the determination of the unknown functions. If one 

attempts to supplement them by obtaining equations from the higher order terms in the original equations, 
more unknowns are also introduced. Therefore we instead restrict our attention to those solutions for which 
pI = 0, i. e., we pressume that in any meteorologically significant solution if p^ is independent of 6 and (j), it 
is also independent of t. Then p^ is determined by the initial data, and by equation 25, so is p^. Thus we 
are left with the four equations 49, 39, 38 and 35 for the four unknown functions v'^,w°, and p^ . 
Equation 35 becomes, since p\ = p\ =0, 

(51) u\pI - ^p°) = 0. 

If the second factor, determined by the initial data, is not zero (implying that the zero order solution is 
non-isentropic) then = 0. Equations 49, 39 and 38 then suffice to determine v'^.uP and p"^. 

On the other hand, if the second factor in equation 51 is zero, implying that the zeroth order solution is 
isentropic, this equation is useless and we remain with three equations for four unknowns. To obtain another 
equation we equate to zero the coefficient of in equation 15 and find 

(52) p u {--—), = -{p^--^p^). 

Now from equations 32 and 33, is independent of 9 and (j>. Therefore we restrict our attention to solutions 
p^ independent of t on the basis of the presumption mentioned above. Then by eq. 31 p^ is also independent 
of t and both p^ and p'^ are determined by the initial data. Equation 52 now becomes 

(53) " " = °- 

Here again the second factor may not be zero, implying the solution is not isentropic to first order, and 
then = 0. Then, as before, equations 49, 39 and 38 suffice for the determination of v^,^]^ and p^. If 
the second factor is zero, equation 53 is useless and we equate to zero the coefficient of in equation 15 to 
obtain the additional equation 

(54) p^ip')--' - -Ip'pW-''-' + »>i(/>°)-' - VII 

2^ IP^ IV^P^ 7(7 + 1) {p^ f ^ 



Now we have the five equations 49, 39, 38, 37, and 54 for the determination of (? and u'^ 
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6 Summary of Results. 



By introducing a certain transformation of variables involving a parameter e, and by assuming that the 
solution can be expanded in powers of e, we have obtained a simplified system of equations for the de- 
termination of the first terms in the expansion of the solution. These simplified equations imply that the 
pressure is hydrostatic and the horizontal wind geostrophic (to the order in e considered). In the course 
of the derivation it was found that and are independent of 6 and 0. We consequently restricted our 
attention to solutions in which these quantities are also independent of t, presuming that any other solutions 
are not of meteorological importance. There are two sets of simplified equations, depending upon the degree 
of isentropy of the initial data. These two sets are considered separately below. 



6.1 Nonisentropic Case. 

This case obtains if at least one of the quantities — and {^ — -^)z is not zero. Then u'^+eu^+e'^u'^ = 

0, p°{z) + ep^{z) + e'^p^{z) is given by the initial data and p^{z) + ep^{z) + e'^p'^{z) + e^p^{9, (/>, z,t) is determined 
by the hydrostatic equation. Equations 38, 39 and 49 determine v'^,w'^ and p'^. These equations involve no 
z derivatives, and if p = p^/p° is introduced as a new unknown, the coefficients are also independent of z. 
The equations then become, omitting the superscript on and w^: 

(55) w = {2pcos6o)~^pg, 

(56) V = {—2ij,cos9osm0o)~^Pfi„ 

(57) {vt + VV0 + wv^[sm.0o]~'^)^[sm0o]~'^ 

-{wt + VW0 + ww^[sm(t)\~'^)e - 2[cos6io]~V<^ = 0. 



6.2 Isentropic Case. 

This case obtains if both — ^p-p" and (^ — -^)z are zero. Then + eu^ = 0, p^{z) + ep^{z) + e^p'^(z) 
is given by the initial data and p'^{z) + €p^{z) + f^p'^{z) + e^p^{9,(p,z,t) is determined by the hydrostatic 
equation. Equations 38, 39, 49 and 54 determine and u^. Omitting superscripts and eliminating 

by means of equation 37, these equations become: 

(58) w; = (2/ip°cos6io)-V, 

(59) w = (—2/1/9° cos 00 sin 0o)~^lJ0j 

(60) 2i^cos0o{p°y\p°u)z + (vt + wo + wv^[sm0or^)4>iin0o]~'^ 

-{wt + vwg + ww^[sin9()]^^)o - 2[p°cos6'o]"V0 = 0, 

^p^ ^P^ _i 

(61) {p + T — ^Pz)t + v{p + - — -pz)e + w[sin6»o] {p + t — ^Pz),j, 

MP Mp MP 

A simplification of these equations results if the coefiicient of u in equation 61 is zero, which may be called 
the extreme isentropic case. Then equations 58, 59 and 61 can be solved for p, w and v and then u can be 
found from equation 60. 



7 Boundary Layer Effect. 

It may be noticed that the initial data must satisfy various conditions, i. e., geostrophic and hydrostatic 
equations. Similar conditions must be satisfied by the boundary data on the spacial boundaries. Stated 
otherwise, all the initial and boundary data cannot be prescribed arbitrarily, as one would have expected. 
This is typical of the boundary layer phenomenon which always arises in the asymptotic expansion of the 
solution of a system of differential equations, because of the reduced order of the approximate system. The 
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question arises as to the; proper choice of data for the approximate solution, when the data for the exact 
problem are given, in order that the approximate solution best approximate the exact solution away from 
the boundaries. This difficult question should not be important in the present case, however, since the 
boundaries are not "real" but are within a larger region in which the asymptotic solution is presumably 
valid. Therefore the initial and boundary data, if obtained from observations, should satisfy the required 
conditions. 



8 The Barotropic Atmosphere. 

If the atmosphere is barotropic, i. e., if there is a functional relation between p and p, then this relation 
replaces the entropy equation, eq, 5 



(5') P = f{p). 

To derive the simplified equations in this case, we proceed exactly as before, but replace all consequences of 
eq. 5 by those of eq. 5'. Thus instead of eqs. 22, 29 and 35 we have 



(22') p° = f{p°), 

(290 p' = rip')p\ 

(35') p2^/Y)[pi]2 + /'(p0y. 

The derivation of eq. 44 is the same as before, but to proceed further we restrict our attention to solutions 
p° independent of t, since is already independent of and by eqs. 19 and 20. Then from eq. 44 we find 
that = 0. We further consider only solutions such that p^ and p^ are independent of t, since both are 
independent of 9 and (f). Equations 38, 39 and 49 follow as before for the determination of w", and u^. 

From equations 38 and 39 we find that v° and are independent of 2;. To show this, we differentiate 
eq. 38 with repect to z: 

(62) 

The last equality follows from eqs. 18 and 37. Now from eq. 5' we have 

(63) PePz - PzP0 = 0. 

Since p",p^ and p^ are independent of 9 and (j), the lowest order term in eq. 63 is p°pg — Pgp^ which is 
consequently zero. Thus from eq. 62, = and similarly v'^ =0. 

Making use of these results, we can eliminate fron eq. 49 by integrating that equation with respect to 
z from to z, obtaining 

(64) pV = + + 



sm Uo sm Uq 

(«;? + v'^wl + - 2p|[p° cos^o]-^} , , . [p\z) - /(O)]. 

SmCg Z/iAiCOSWo 

In eq. 64 we have made use of eqs. 18 and 36. Now aX z = oo the left side vanishes, and since p°((X)) also 
vanishes while p°(0) is positive, the expression in braces on the right must vanish. This is just equation 57. 

Since this expression is independent of z, we find from eq. 64 that = 0. Therefore the quantities tP,uP 
and p'^ are determined by eqs. 55, 56 and 57 in the barotropic atmosphere, just as in the non-isentropic case 
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(Subsection 6.1) for a baroclinic atmosphere. This residt is somewhat surprising, since one might expect 
the barotropic atmosphere to correspond to the isentropic case. This is also the case, however, for if one 
assumes that the atmosphere is exactly isentropic (or at least is up to third order) then eqs. 58 and 59 imply 
that and uP are independent of z; eq. 60 then yields eq. 57 and the result = as above, and eq. 61 
becomes an identity. 

9 Special Solutions (Nonisentropic Case). 

9.1 Zonal Motion: v — — 

Equations 55 and 57 yield 

Peet = 0. 

Thus, with a, b, c arbitrary functions, we have 

p = a{t)0 + b{t) + c(9), 

w = {2ii cos 0o)-^ [a{t) + c' {9)]. 

9.2 Meridional Motion: w — pe — 

Equations 56 and 57 yield 

Pcjxpt + 4/U sin^ doP4> = 0. 

Integrating 

p^t + 4/xsin^ 9qP = a{t). 

Thus we have 

p = ai{t)+ ( /(a)e"'^-(^''^'"'^°/")*da, 
Jc 

V = (-2/icos^osin0o)"'y" Q:/(a)e"'^-(^''"'"'^°/"^*rfa. 
If we impose periodicity in cj), then a = n (n = 0, ±1, ±2, • • •) and the integral is replaced by a series. 

9.3 Perturbation of Zonal Motion. 

Assume a solution analytic in a parameter r] which yields a steady zonal motion for 77 = 0. The solution 
may be written: 

p = p''{0) + vp\e,c^,t) + 7j^p^ + ---, 

w = w"{9) +r]pl{2ficoseo)~'^ -\ , 

V = (-^ sin 26*0 sin 00)^^00* + w" (6*) pj,^^(-^ sin 26*0 sin^ 6*0)"^ - (2/xcos6'o)"VJ,0t 
+ inam29o)-\w°pl)e - (2/xcos0o sin9o)-\w°pl^)e - 2{cos9o)-'pl = 0. 

li = this simplifies still further to 

-r^pl4,4, + Pm + 4mpJ, = 0. 
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